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Abstract 

We present an analysis of the classic Alcubierre metric based on conformal gravity, rather than 
standard general relativity. The main characteristics of the resulting warp drive remain the same as 
in the original study by Alcubierre, namely that effective super-luminal motion is a viable outcome 
of the metric. 

We show that for particular choices of the shaping function, the Alcubierre metric in the context 

cr 

of conformal gravity does not violate the weak energy condition, as was the case of the original 
solution. In particular, the resulting warp drive does not require the use of exotic matter. 

Therefore, if conformal gravity is a correct extension of general relativity, super-luminal motion 
via an Alcubierre metric might be a realistic solution, thus allowing faster-than-light interstellar 
travel. 
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I. INTRODUCTION 

In 1994 M. Alcubierre introduced the so-called Warp Drive Metric (WDM), within the 
framework of General Relativity (GR), which allows in principle for super- luminal motion, 
i.e., faster-than-light travel [l|. This super-luminal propulsion is achieved by respectively 
expanding and contracting the space-time behind and in front of a spaceship, while the 
spacecraft is left inside a locally flat region of space-time, within the so-called warp bubble. 

In this way the spaceship can travel at arbitrarily high speeds, without violating the laws 
of special and general relativity, or other known physical laws. Furthermore, the spacecraft 
and its occupants would also be at rest in flat space-time, thus immune from high acceler- 
ations and unaffected by special relativistic effects, such as time dilation. Enormous tidal 
forces would only be present near the edge of the warp bubble, which can be made large 
enough to accommodate the volume occupied by the ship. 

However, Alcubierre [l] was also the first to point out that this hypothetical solution 
of Einstein's equations of GR would violate all three standard energy conditions (weak, 
dominant and strong; see j^j, {3)], and ^| for definitions). In particular, the violation of the 
weak energy condition (WEC) implies that negative energy density is required to establish 
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the Alcubierre WDM, thus making it practically impossible to achieve this type of super- 
luminal motion, unless large quantities of exotic matter (i.e., with negative energy density) 
can be created. Since our current knowledge of this type of exotic matter is limited to some 
special effects in quantum field theory (such as the Casimir effect), it is unlikely that the 
Alcubierre WDM can be practically established within the framework of General Relativity. 



) or refining and 



Following Alcubierre's seminal paper, many other studies appeared in the literature, 
either proposing alternatives to the original warp drive mechanism ( p | , [6 | ) 
analyzing in more detail the original idea (Q, {9], 10], 11], [l2|, [l3|, Q, [15|, [16|, 
\v\ . jl8|. 19]). However, all these studies were conducted using standard GR and could not 



avoid the violation of the WEC, meaning that some exotic matter would always be required 
for faster-than-light travel. Similar issues also exist in other well-known GR solutions for 
super-luminal motion, such as space-time wormholes (20 ]. 

Einstein's General Relativity and the related "Standard Model" of Cosmology have been 
highly successful in describing our Universe, from the Solar System up to the largest cos- 
mological scales, but recently these theories have also led to a profound crisis in our under- 
standing of its ultimate composition. From the original discovery of the expansion of the 
Universe, which resulted in standard Big Bang Cosmology, scientists have progressed a long 
way towards our current picture, in which the contents of the Universe are today described 
in terms of two main components, dark matter (DM) and dark energy (DE), accounting for 
most of the observed Universe, with ordinary matter just playing a minor role. 

Since there is no evidence available yet as to the real nature of dark matter and dark 
energy, alternative gravitational and cosmological theories are being developed, in addition 
to standard explanations of dark matter/dark energy invoking the existence of exotic new 
particles also yet to be discovered. In line with these possible new theoretical ideas, Con- 
formal Gravity (CG) has emerged as a non-standard extension of Einstein's GR, based on 
a possible symmetry of the Universe: the conformal symmetry, i.e., the invariability of the 
space-time fabric under local "stretching" of the metric (for reviews see 21], j^j). This 
alternative theory has been re-introduced in recent years (following the original work by H. 



Weyl 



23]. JjSl 



DM and DE 



eading to cosmo 



26|, |27j, |28|, |29j, |30j, 



ogical models which do not require the existence of 



3l|, |32j). 



In view of a possible extension of Einstein's General Relativity into Conformal Gravity, 
in this paper we have re-considered the Alcubierre WDM, basing it on CG rather than 
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standard GR. In Sect. [IT] we review the fundamental principles of CG and the calculation of 
the stress-energy tensor in this gravitational theory. In Sect. II I II we consider the Alcubierre 
metric in CG and compute the energy density for different shaping functions of the metric. 

In particular, we will show that for certain shaping functions, the Alcubierre metric in 
the context of Conformal Gravity does not violate the weak energy condition, as was the 
case of the original solution. This analysis continues in Sect. IIV} where we study other 
energy conditions and estimate the total energy required for this CG warp drive. Finally, in 
Sect. |Vl we conclude that if CG is a correct extension of GR, super-luminal motion via an 
Alcubierre metric might be a realistic possibility, thus enabling faster-than-light interstellar 
travel without requiring exotic matter. 



II. CONFORMAL GRAVITY AND THE STRESS-ENERGY TENSOR 



H. Weyl in 1918 ( 23|], 24|, 25|) developed the "conformal" generalization of Einstein's rel- 
ativity by introducing the conformal (or Weyl) tensor, a special combination of the Riemann 
;ensor Rx^k, the Ricci tensor R^ u = R x ^\ v and the curvature (or Ricci) scalar R = 



33): 



C\^ UK = Rxhuk — -^(g\vRnK — 9\kR^v — g^R\K + g^K,R\u) + ^R {g^g^n — gx^^u), (l) 

where C x ^\ u {x) is invariant under the local transformation of the metric: 



9»u( x ) ->• 9^u(x) = e 2a(x) gflu {x) = Q 2 {x) gflu {x). (2) 

The factor = e a<yX ^ determines the amount of local "stretching" of the geometry, hence 
the name "conformal" for a theory invariant under all local stretchings of the space-time 
(see 22( and references therein for more details). 

This conformally invariant generalization of GR was found to be a fourth-order theory, as 
opposed to the standard second-order General Relativity, since the field equations originating 
from a conformally invariant Lagrangian contain derivatives up to the fourth order of the 



metric, with respect to the space-time coordinates. Following work done by R. Bach 34j . 



C. Lanczos 



35[ and others, CG was ultimately based on the Weyl or conformal action: 1 



1 In this paper we adopt a metric signature (-,+,+,+) and we follow the sign conventions of Weinberg [3 



' w 



-OLc 



d 4 x {-g) 1 ' 2 C\fjiUK C X ^ K , 



(3) 



or on the following equivalent expression, differing from the previous one only by a topolog- 
ical invariant: 

(4) 



Iw = ~2a g J d 4 x (-g) 1 ' 2 (r^R^ - ^R 2 ^j 



where q = det(q,,, v ) and a g is the gravitational coupling constant of Conformal Gravity (see 



39j). 2 Under the conformal transformation in Eq. ([2]), the Weyl tensor 
transforms as C\^ VK — > C\^ VK = e 2a<yX ^C\^ vli = £l 2 (x)C\^ UK , while the conformal action Iw is 
locally conformally invariant, the only general coordinate scalar action with such properties. 
R. Bach [34J introduced the gravitational field equations in the presence of a stress-energy 



tensor 3 T, 



W, 



T 



(5) 



as opposed to Einstein's standard equations, 



G 



(J.V 



R/J.U c^Qiiv R 



8ttG 



T, 



(6) 



where the "Bach tensor" W^ u [34] is the equivalent in CG of the Einstein curvature tensor 
G^y on the left-hand side of Eq. (J6]). 



Wn U has a very complex structure and can be defined in a compact way as 



4Q|: 



— 2C a ^-^. a + C° ' ^ R/3a, 



(7) 



or in an expanded form as ([38], |4l|): 



In this section we will leave fundamental constants, such as c and G, in all equations, but later we will 
use geometrized units (c = 1, G = 1), or c.g.s. units when needed. 

2 In these cited papers, a g is considered a dimensionless constant by using natural units. Working with 
c.g.s. units, we can assign dimensions of an action to the constant a g so that the dimensionality of Eq. 
([3]) will be correct. 

3 We follow here the convention [36( of introducing the stress-energy tensor T^ v so that the quantity cXbo 
has the dimensions of an energy density. 
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~Wflv = —-g^iv R ,X ;X + -^R;tx,P + R/mv' ;\ ~ Rp ;v;\ ~ Rv -,n;\ + ~^R RfjLV (8) 

— 2R fl x R\ u + -g^v R\ p R xp — -gnu R 2 , 

involving derivatives up to the fourth order of the metric with respect to space-time coordi- 
nates. 

Therefore, in Conformal Gravity, the stress-energy tensor is computed by combining 
together Eqs. ()5]) and ([8]): 



1 2 

= ^OL g W^y = Aa g (—-g^u R' ; a + -^R-^u + Rpu ; a — Rp -,v,\ — Ru -^-x (9) 

2 1 i 

+ -^R Rfj,v — 2i? M A R Xv + -g^ u R\ p R xp — -g^ R 2 ). 

This form of the tensor will be used in the following sections, in connection to the Alcubierre 
metric, to compute the energy density and other relevant quantities. 

For this purpose, we have developed a special Mathematica program which enables us to 
compute all the tensor quantities of both GR and CG, for any given metric. In particular, 
this program can compute the conformal tensor Ca mj , k in Eq. (TJJ), the Bach tensor in 
Eq. (jSJ), or the stress-energy tensor T^ v in Eq. ([9]) by performing all the necessary covariant 
derivatives. Given the complexity of these types of computations, we have tested extensively 
our program against the results for C\ llVK and W^ v computed by Mannheim et al. in Ref. 



391 ] for different metrics, obtaining a perfect agreement. 



III. ALCUBIERRE METRIC, SHAPING FUNCTIONS AND THE WEAK EN- 
ERGY CONDITION 

The original Alcubierre metric jjj considered a spaceship traveling along the x-axis, with 
motion described by a function x s (t) and spaceship velocity v s (t) = dx ^ . Using the 3+1 
formalism of GR, the metric was written in Cartesian coordinates as (c = 1): 

ds 2 = -dt 2 + [dx - v s (t) f{r s ) dt] 2 + dy 2 + dz 2 , (10) 
where r s is the distance from the spaceship position: 



6 



r,(t) = ^[x-x s (t)} 2 + y 2 + z 2 , (11) 

and f(r s ) is a "form function" or "shaping function" which needs to have values / = 1 and 
/ = respectively inside and outside the warp bubble, while it can have an arbitrary shape 
in the transition region of the warp bubble itself. 

The original shaping function used by Alcubierre was: 

_ tanh [a (r s + R)} - tanh [a (r s - R)} 
I[rs) ~ 2tanh(a J R) ' 1 ] 

where R > basically indicates the radius of the spherical warp bubble, while a > relates 
to the bubble thickness, which decreases with increasing values of a. In the following, we 
will refer to the function in Eq. (112]) as the "Alcubierre shaping function" (ASF). 

We will show that the particular form of the shaping function can play an important role 
in the energy conditions for the WDM. In our analysis we tested several different functions 
obeying the general requirements for / outlined above. In addition to the Alcubierre function 
above, in this paper we will also use the following: 

m = < - ( » r ; r ' < i , d3) 

; r s > Rj 

where m is a positive integer. Since this particular function for m = 4 is used by J. Hartle 



to illustrate the warp drive in his textbook 20j, we will refer to the function in Eq. ({TBI) as 
the "Hartle shaping function" (HSF). 

The top panels in Fig. 1 illustrate the differences between the Alcubierre shaping function 
(top left panel, for a = 8 as used originally by Alcubierre) and the Hartle shaping function 
(top right panel, for m — 4). All functions in this figure are computed for a fixed value of 
R = 1 and at time t — 0, when the spaceship is located at the origin. All quantities shown 
in the different panels are plotted as a function of the x coordinate of the spaceship motion 



and of the transverse cylindrical coordinate p = ^/y 2 + z 2 . Similar coordinates will also be 
used in the other figures. 4 



4 The cylindrical coordinate p = \Jy 2 + z 2 should be considered as non-negative and all quantities in the 
figures plotted only for p > 0. However . for illustrative purposes and also to follow similar figures in the 
literature (such as those in 12|, etc.), we decided to let p run on negative values in all figures, 

except in the last one, where we restrict p > for a correct energy calculation. 
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The expansion/contraction function 9 of the volume elements behind/in front of the 
spaceship was also computed by Alcubierre as 

9 = v s (t) [X ~ Xsit)]df } rS \ (14) 

and is illustrated for v s = 1 (in geometrized units, i.e., v s = c in traditional units) in the 
second row of Fig. 1, for the two different shaping functions. Again, the choice of the 
parameters R, a, and m is the same as in the top panels in the figure. The expansion 9 for 
the ASF is the same as Fig. 1 in Ref. 1], while the corresponding 9 for the HSF is slightly 
different, but still shows expansion of the normal volume elements behind the spaceship and 
contraction in front of it. 

The weak energy condition ) requires that T^tH" > for all timelike vectors 

t M . Alcubierre has also shown that for the Eulerian observers in the warp drive metric, and 

n 

for their 4- velocity n M , the following relation holds [1J: 

,00 1 V sP 2 ( d f 



s 

->oo 



Vn^ = r~ = --^ , (15) 



which implies that the energy density T is negative everywhere for any choice of the shaping 
function / and, therefore, the WEC is violated (also the dominant energy condition -DEC- 
and the strong energy condition -SEC- are violated in the analysis based on GR 

This violation of the WEC in GR is illustrated in the third row of Fig. 1, where T 00 
is calculated using Eq. (JTBI) for both shaping functions. Although the results in the two 
panels are slightly different, they obviously show negative energy densities and therefore a 
complete violation of the WEC. 

The situation is different if we compute the energy density T 00 in the framework of CG, 
following Eq. fl9]), setting a g — 1 for simplicity, and using the completely contravariant 
form of the stress-energy tensor, instead of the covariant one. As seen in the bottom row of 
Fig. 1, the energy density in CG is completely different from the one calculated within the 
framework of GR. In the bottom left panel T 00 is computed with the ASF and the resulting 
function is in part positive and in part negative, thus still violating the WEC. 

However, the bottom right panel shows T 00 computed with the HSF and in this case the 
energy density is completely non-negative, showing that the WEC is verified and no exotic 
matter is needed to establish the warp drive. This non-negative energy density plot in the 
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bottom right panel of Fig. 1 is the main result of our paper as it shows that - if CG is the 
correct extension of GR - it might be possible to establish a warp drive without having to use 
negative energy (mass), thus overcoming the main difficulty of the warp drive mechanism. 

The explicit expression of T 00 in CG, computed with our Mathematica program, is rather 
cumbersome and is reproduced in Eq. (I A . 1 1) of the Appendix. Here we present just the 
graphical computation of T 00 in Fig. 1, or in the other figures in this paper. We have 
tested the validity of these results by running the program in several different ways, includ- 
ing computing the stress-energy tensor for a simplified three-dimensional Alcubierre metric 
(coordinates x, y, t only), always obtaining consistent results. The differences in the two 
plots at the bottom of the figure can be attributed to the different shaping functions and 
their derivatives up to the fourth order. All these derivatives enter the complex expression 
of T^ v in CG, as in the master Eq. fl9]), and their interplay ultimately determines the shape 
of T 00 , or of the other components, in a way which is hard to predict before the actual 
computation is performed. 

In figure 2 we analyze the dependence of the energy density T 00 on the spaceship velocity 
v s . In this case we consider only the Hartle shaping function for m = 4 and R = 1, and 
we compute T 00 in CG (a g = 1) for speeds ranging from the sub-luminal v s = 0.25 c to the 
super- luminal v s = 3.00 c (for v s = 1.00 c we obtain the same function as in the bottom 
right panel of Fig. 1). The shape of the energy density function is about the same for speeds 
up to v s = 1.50 c, although the function values increase with speed. For higher velocities, 
the function develops two "downward lobes" which eventually become negative for speeds 
v s > 2.50 c. This implies that the WEC is verified for speeds up to v s ~ 2.50 c, while at 
higher velocities exotic matter would be required to sustain the warp drive. 

This apparent "speed limit" at about v s ~ 2.50 c might be raised or overcome completely 
by adopting a different shaping function, instead of the HSF used here, but this analysis 
would go beyond the scope of this work. In any case, the results reported in Fig. 2 show that 
a warp drive in CG with positive energy density is possible for a wide range of spaceship 
velocities; therefore, if CG is the correct extension of GR, the Alcubierre warp drive might 
be a viable mechanism for super-luminal travel. 

In Fig. 3 we present the other components of the stress-energy tensor. These were 
computed with the same Mathematica program, following Eq. fl9]) with a g — 1, leading to 
even more complex expressions than the one for T 00 (we will omit to report these expressions 
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for brevity). To simplify the computation, we used the covariant components T^ u and 
adopted cylindrical coordinates around the x-axis, (t,x,p,<p) = (0,1,2,3), instead of the 
Cartesian coordinates of the original Alcubierre metric. The results in the different panels 
are labeled accordingly. We recall that the stress-energy tensor is symmetric, T^ v = T Ufl , and 
only the non-zero components are illustrated in this figure, under similar conditions used 
before (CG with HSF and parameters m = 4, R = 1, v s = 1.00 c). 

The shapes of the other components of T^ v are more complex than the one of T 00 , but 
they can all be determined analytically, either in covariant or contravariant form, using our 
program. If this exact form of the stress-energy tensor could be established in the region 
surrounding the spacecraft, warp drive motion would be feasible within the framework of 
Conformal Gravity. 

We also want to point out that we set the spaceship motion in the positive direction of 
the x-axis (setting v s = +1.00 c), and this is reflected in the symmetry, or lack thereof, of the 
components of T^ u . While some components, such as T 00 (or Too), T i = T w , Tn, T 2 2, T 33 , 
appear to be symmetric under the exchange x — > —x, the other components, Tq 2 = T20 and 
T12 = T21, are not symmetric under this exchange and, therefore, these components must 
contain information about the spaceship direction of motion. This argument overcomes the 
objection addressed in Ref. [3] that since T 00 is symmetric about the x s = plane, there 
is uncertainty in where the space-time is expanded/contracted, thus making it impossible 
for the spaceship to know in which direction of the x-axis, positive or negative, to move. 
Rather, Fig. 3 shows that a "bias" towards one of the two possible directions is induced by 
some components of T^ u . 

Fig. 4 illustrates one last dependence of the energy density T 00 on the parameters used. 
In this case, we set R — 1, v s — 1.00 c and consider the Hartle shaping function as in Eq. 
(fl3|) . while varying the integer parameter m. In addition to our standard value, m = 4, we 
have also tried values from m = 2 to m = 10, as shown in the figure. 

In general, increasing the m value increases the internal volume of the warp bubble, 
where space-time is flat, and therefore decreases the thickness of the bubble wall where the 
space-time distortion takes place, in a way similar to that of the a parameter in the original 
Alcubierre warp drive. However, increasing m also increases the energy required to establish 
the warp drive (as we checked by integrating the functions in Fig. 4). Therefore, it appears 
to be convenient to use a low value for this parameter. As shown in the different panels 
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of the figure, the first value, m = 2, does not work since it does not leave a flat space- 
time volume inside the bubble, while m = 3 seems to create a very small volume inside the 
bubble. Therefore, intermediate values such as m = 4 — 6 would appear to be more adequate 
to establish the warp drive. Meanwhile, the solution for m — 10 would require much more 
energy and would not give any advantage, except reducing the thickness of the warp bubble. 
We also checked that changing the value of m does not have a strong effect on the "speed 
limit" of v s ~ 2.50 c, reported above for the case m — 4. Thus, this value of the parameter 
seems to be the most adequate for this type of solutions. 

IV. OTHER ENERGY CONDITIONS AND WARP DRIVE ENERGY ESTIMATE 

In the previous section we have discussed at length the weak energy condition -WEC- 
for the Conformal Gravity Alcubierre warp drive. We have seen that, if the Hartle shaping 
function is used, this condition is not violated for a wide range of spaceship velocities, 
including super-luminal speeds. In this section we will briefly analyze the other main energy 
conditions and estimate the energy necessary to establish the warp drive in CG. 

The dominant energy condition -DEC- is reported in the literature 
\T^\ for any //, u, or equivalently as assuming the WEC plus the additional condition that 
T^t^ is a non-spacelike vector, i.e., T^ U T U < 0. It is easy to see that using as a vector 
the 4- velocity of the Eulerian observers [1], the previous condition for the DEC becomes 
T° A T 0A < 0. 

Figure 5 illustrates the violation of the DEC for our standard solution (AWD with HSF 
and m = 4, R = 1, v s = 1.00 c). The plotted function T° aT 0A is not negative everywhere, as 
required by the DEC, but shows a violation for the central portion of the warp bubble. Even 
if this energy condition appears to be violated, this does not notably affect the feasibility of 
our CG warp drive. We recall that the DEC is usually related to the standard perfect fluid 
stress-energy tensor, T^ u — (p + p)U^U u + pg^, where here p and p are the fluid density 
and pressure, while is the fluid 4- velocity. In this context the DEC requires p > \p\, but 
this condition is not required in general by all classical forms of matter therefore, its 
violation in our case is not particularly significant. 

On the contrary, our standard solution also verifies the strong energy condition -SEC- 
which states that T^tH" > \T X x t a t a for all timelike vectors V. A gain, using the Eulerian 
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4- velocity vector in place of the previous condition is equivalent to T 00 + |T A \ > 0. Since 
the scalar T x \ is identically zero for all our solutions, as checked using our Mathematica 
program, the SEC is equivalent to T 00 > 0, which is the WEC already verified in Sect. Ill IL 
Finally, we want to estimate the energy necessary to establish our CG warp drive, under 
reasonable conditions. For this purpose, in Fig. 6 we computed once again the energy 
density T 00 for our AWD with the Hartle shaping function (a g = 1, m = 4, v s = 1.00 c), 
but this time for R = 10000 cm = 100 m, a reasonable radius for a warp bubble enclosing 
our spaceship. 

Figure 6 illustrates this solution, plotted only for p = \J y 2 + z 2 > 0, as this is the correct 
interval for the transverse coordinate p. The cylindrical symmetry of this solution can also 
be better appreciated in this type of plot. We then followed the procedure outlined in 



Refs. 



,8 



llj . to integrate the local energy density over the proper volume, in cylindrical 



coordinates at time t = over all space, obtaining the total energy E: 

E = a g c J d 3 xy/\g\ T 00 = (1.86 x 10 10 a g ) erg, (16) 

where g = Det \gtj\ is the determinant of the spatial metric on the constant time hyper- 
surface. Since we assume that the spaceship is traveling at constant velocity, v s = 1.00 c, 
the total energy is also constant with time. In the last equation, we reinstated a factor of c 
to obtain the correct dimensions (see footnote before Eq. (JSJ)) and also inserted an overall 
multiplicative factor a g , which corresponds to the conformal gravity coupling constant in 
Eq. ()9]). This factor is necessary since our computation of T 00 in Fig. 6 was done assuming 
a g = 1. 

Therefore, we need to know the CG value for a g in order to complete our energy esti- 
mation. Unfortunately, the value of this coupling constant is not well determined yet. The 



only value in the literature is reported by P. Mannheim 



30|: 



a g = 3.29 x 10 82 ergs, (17) 

since this coupling constant has the dimensions of action. Inserting this value for a g in Eq. 
f fT6|) . we obtain the energy estimate: 

E = 6.12 x 10 92 erg = 6.81 x 10 71 g = 3.42 x 1O 38 M , (18) 
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where we also converted this energy into equivalent mass and compared our result with the 
Solar mass M & = 1.99 x 10 33 g. 

The estimate in Eq. (I18p would imply that an enormous amount of (standard) mass- 
energy is needed to establish our warp drive at a velocity equal to the speed of light, with 
a reasonable size for the warp bubble. However, the CG value of a g is not well-established, 
since the number in Eq. (1171) represents only an estimate of the macroscopic value of this 
coupling constant. This does not need to be the same as the microscopic value associated 
with the fundamental theory, which could be reduced by a factor of N, where iV could be 
the number of occupied baryonic states in a galaxy (N ~ 10 68 ), or possibly the number of 
baryons in the Universe (N ~ 10 80 ) 42]. 

Therefore, our estimate could be reduced by many orders of magnitude. Moreover, the 
energy necessary to establish the warp drive might also be decreased by using a more efficient 
shaping function, an analysis which we leave for a future study on the subject. 



V. CONCLUSIONS 



In this paper we have analyzed in detail the Alcubierre warp drive mechanism within 
the framework of Conformal Gravity. We have seen that a particular choice of the shaping 
function (Hartle shaping function, instead of the original Alcubierre one) can overcome the 
main limitation of the AWD in standard General Relativity, namely the violation of the 
weak energy condition. 

In fact, we have shown that for a wide range of spaceship velocities, the CG solutions do 
not violate the WEC, and, therefore, the AWD mechanism might be viable, if CG is the cor- 
rect extension of the current gravitational theories. All the components of the stress-energy 
tensor can be analytically calculated, using a Mathematica program based on Conformal 
Gravity. Thus, a warp drive can, at least in principle, be fully established following our 
computations. 

We have also checked two other main energy conditions: the SEC is always verified, while 
the DEC is violated, at least in the case we considered. Finally, we estimated the energy 
needed to establish a reasonable warp drive at the speed of light. This energy depends 
critically on the value of a g , the CG coupling constant, which is not well known. Therefore, 
this estimate will need to be refined in future studies. 
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VI. APPENDIX: ENERGY DENSITY EXPRESSION IN CONFORMAL GRAV- 
ITY 

We present here the expression for the energy density T 00 in Conformal Gravity, com- 
puted using our Mathematica program. This is the general form of T 00 for any shaping 
function f[r s ] and its derivatives, up to the fourth order. The energy density is a function 
of coordinates t, x, p = yfy 2 + z 2 , and, therefore, it has a cylindrical symmetry around the 
x-axis. The distance r s is defined in Eq. ( II ip . v s is the spaceship velocity and a g is the 
Conformal Gravity coupling constant. 



T 



oo 



'3tv, 



Ar s v s 2 p 2 (Q(x - v s t) 2 + 5p 2 ) (-1 + f[r s })f[r 



(A.l) 



+ Ar A s v 2 s ((x - v s t) 2 + 3p 2 ) f[r s ] A + f'[r s } 2 (-24(x - v s tf + 3 (l + 5v s 2 ) (x - v s t) 2 p 2 

+ (27 + 10v s 2 ) p 4 + v s 2 (5p 2 (3(x - v s t) 2 + 2p 2 ) (-2 + f[r s \)f[r s \ - Ar 2 {x - v s t) 2 

(A(x - v s tf + 3p 2 ) (-1 + f[r s \)f"[r s \)) - 2r s f'[r s }(lQ(x - v s tf - 8p 2 + (-60 (a; - v s tf 

+ (-97 + 5v 2 ) {x - v s t) 2 p 2 + (-37 + 3^ s 2 ) p 4 )f"[r s } + r s (4 (-4 + v 2 ) (x - v s t) 4 

+ 3 (-9 + v 2 ) (x - v s tfp 2 - llp 4 )/^[r s ] + v 2 f[r s ] 2 (p 2 (5(z - ^t) 2 + 3p 2 ) f [r s ] 

+ r s {x - v s t) 2 (4(x - v s t) 2 + 3p 2 ) / (3) [r s ]) + 2/[r s ](-8(x - v s t) 2 + 4p 2 + v s 2 

(-P 2 (5(a - v s t) 2 + 3p 2 ) /" [r s ] - r s {x - v s t) 2 (A{x - v s t) 2 + 3p 2 ) / (3) [r s ]))) 

+ r s 2 (-16 (2(x - v s t) 2 - p 2 ) (-1 + f[r s ])f"[r s ] + (4 (6 + ^ s 2 ) (x - v s t) 4 + (43 + 3^ s 2 ) 

(x - v s t) 2 p 2 + 19p 4 + v s 2 (x - v s t) 2 {A(x - v s t) 2 + 3p 2 ) (-2 + f[r,])f[r,])f" [r s ] 2 

+ 8r.(-l + /[r.]) ((2(x - ^t) 2 + p 2 ) /( 3 )[r s ] + r s (x - v.f) 2 / (4) [r.]))) 

In the main part of our work, we used the Hartle shaping function in Eq. (|13|) . or more 
explicitly: 
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where [...] indicates the positive part of the function. The derivatives of the HSF, up to 
the fourth order, were computed in terms of the Heaviside step function H(x) and the Dirac 
delta function 5(x), also using the following relations for the derivatives of these special 
functions: ^ = sgn(x) = 2H(x) - 1; ™ = 5(x); x n ^^ = (-l) n n\ 5{x). 
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By inserting these derivatives into Eq. ( 1A.1I) . we obtained the expression used to compute 
T 00 in Fig. 1 (bottom right panel), Fig. 2, Fig. 4, and Fig. 6. 
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FIGURE CAPTIONS 



Figure 1: Results for the two different shaping functions (left column ASF, right column 
HSF), computed with parameters v s — 1, R — 1, a — 8, m — 4, a g — 1, t — 0. Top row: 
Alcubierre and Hartle shaping functions. Second row: expansion of the volume elements 
6. Third row: energy density T 00 computed with General Relativity. Bottom row: energy 
density T 00 computed with Conformal Gravity. The WEC is verified in the case shown in 
the bottom right panel. 

Figure 2: Energy density T 00 computed with Conformal Gravity and the Hartle shaping 
function. Parameters used: R = 1, m — 4, a g — 1, t — 0, and variable v s = 0.25 c — 3.00 c. 
The energy density becomes in part negative for v s > 2.50 c, so the WEC is verified for 
speeds up to about v s ~ 2.50 c. 

Figure 3: Stress-energy tensor components T^, in cylindrical coordinates (t,x,p,(j>) = 
(0, 1,2,3), computed with Conformal Gravity and the Hartle shaping function. Parameters 
used: v s — 1, R — 1, m — 4, a g — 1, t — 0. 

Figure 4: Energy density T 00 computed with Conformal Gravity and the Hartle shaping 
function. Parameters used: v s — 1, R — 1, a g — 1, t — 0, and variable m = 2 — 10. In all 
cases the WEC is verified. 

Figure 5: Violation of the DEC in the case analyzed (CG with HSF and m — 4, R — 1, 
v s = 1.00 c, a g — 1, t — 0). The plotted function T° X T 0X is not negative everywhere, as 
required by the DEC, but shows a violation for the central portion of the warp bubble. 

Figure 6: Energy density T 00 computed with Conformal Gravity and the Hartle shaping 
function, plotted for p > 0. Parameters used: v s — 1, R — 10000, m — 4, a g — 1, t — 0. 
Integrating this local energy density over all space, we obtain an estimate for the total energy 
E required to establish the warp drive. 
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